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Abstract While the Helioseismic and Magnetic Imager (HMI) onboard the 
Solar Dynamics Observatory (SDO) provides Doppler velocity [V], continuum 
intensity [Ic], an d line-depth [Ld] observations, each of which is sensitive to 
the five-minute acoustic spectrum, the Atmospheric Imaging Array (AIA) also 
observes at wavelengths - specifically the 1600 and 1700 Angstrom bands - 
that are partly formed in the upper photosphere and have good sensitivity to 
acoustic modes. In this article we consider the characteristics of the spatio- 
temporal Fourier spectra in AIA and HMI observables for a 15-degree region 
around NOAA Active Region 11072. We map the spatio-temporal-power distri- 
bution for the different observables and the HMI Line Core [II], or Continuum 
minus Line Depth, and the phase and coherence functions for selected observable 
pairs, as a function of position and frequency. Five- minute oscillation power in 
all observables is suppressed in the sunspot and also in plage areas. Above the 
acoustic cut-off frequency, the behaviour is more complicated: power in HMI Iq 
is still suppressed in the presence of surface magnetic fields, while power in HMI 
II and the AIA bands is suppressed in areas of surface field but enhanced in an 
extended area around the active region, and power in HMI V is enhanced in a 
narrow zone around strong-field concentrations and suppressed in a wider sur- 
rounding area. The relative phase of the observables, and their cross-coherence 
functions, are also altered around the active region. These effects may help us 
to understand the interaction of waves and magnetic fields in the different layers 
of the photosphere, and will need to be taken into account in mult i- wavelength 
local helioseismic analysis of active regions. 
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1. Introduction 

Helioseismology relies on the observation of waves, as they affect light from the 
outer layers of the Sun, to infer the structure and dynamics of the otherwise 
invisible deeper layers. By studying the oscillations in multiple wavelengths of 
light - and hence at different heights in the atmosphere - simultaneously, it 
is possible to probe the layers of the atmosphere. Such investigations, aiming 
both to better understand the behaviour of the waves and their interaction with 
magnetic fields, and to improve the inference of subsurface properties, have a 
long history. 

The launch in February 2010 of the Solar Dynamics Observatory (SDO), 
carrying both the Helioseismic and Magnetic Imager (HMI) and the Atmo- 
spheric Imaging Assembly (AIA), provides new opportunities for cross-spectral 
helioseismic analysis, with full-disc, high-cadence images at many UV and EUV 
wavelengths from AIA as well as photospheri c Doppler- ve l ocity, continuum, and 
magnetic data from HMI. As discussed by iHowe et all ([201 lh . the AIA 1600 



A and 1700 A near-ultraviolet bands show a clear signature of the five-minute 
acoustic spectrum that - at least in the Sun-as-a-star case - is much less con- 
taminated by granulation "noise" than the continuum intensity in HMI's visible 
6173 A line. 

Both the effect of surface magnetic activity on helioseismic waves in local 
areas and the phase and coherence relationships between the velocity and the 
intensity of radiation at different wavelengths have been studied for nearly two 
de cades. 

Brown et all ( 1992h reported finding small areas associated with active regions 



that produced a disproportionate amount of acoustic power in the 5.5-7.5 mHz 
frequency band in ground-bas ed Doppler observations using the Fe I 5576 A 
line, while iBraun et all ( 1992 ) made the first observations of defined acoustic 



haloes around active regions in the same frequency band using Ca K inten- 
sity. However, atmosp heric seeing can have confusing effects on ground-based 
intensity observations (Hil l et al , l200lf h particularly in sunspot umbrae. The 



picture becam e clear er with the launch of the Michelson Doppler Imager (MDI: 



Scherrer et a/.LIl995f ) onboard the Solar and Uelio spheric Observatory (SOHO), 
which allowed full-disc o bservations from space in both intensity and velocity 
using the 6768 A Ni line. Ijain and Haberl ((2002) studied the power distribution 
in the line-depth, Doppler-velocity, and continuum-intensity observables from 
MDI. In both sunspot and plage regions they found that the power in areas of 
strong magnetic field was suppressed in all observables, while for velocity and 
line depth (but not for continuum intensity) there was a "halo" of enhanced 
power surrounding the magnet ic-field concentrations, which they concluded was 
acoustic in origin. The Transition Region and Coronal Explorer (TRACE) made 
possible high-resolution, spac e -base d observations in the 1600 and 1700 A bands. 
Judge, Tarbel l, and Wilhelml (|200lf ) used data from another SOHO instrument, 



Solar Ultraviolet Measurements of Emitted Radiation (SUMER), TRACE, and 
MDI to investigate the oscillations in the chromosphere. This work established 
that the chro mospheric modes w ere primarily the same p-modes seen in the 



photosphere. iKrijger et all ([20011 ) used TRACE to study oscillations in the 1600 
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Figure 1. Spatial variation of each observable averaged over 23 May 2010: HMI V (a), HMI 
I c (b), AIA 1600 A (c), AIA 1700 A (d), HMI I L (e), and HMI longitudinal magnetic-field 
strength (f). 



A UV band in the quiet Sun and found enhanced power at three- minute periods 
around patches of network field. iMuglachl ([2003) also looked at acoustic power 
around an active region in short (two - four hour) stretches of TRACE data and 
did not find a high-frequency acoustic halo in the 1600 or 1700 A bands. 

The Hinode/ Solar Optical Telescope has also been used for multi- wavelength 
seismology using the Ca H line and the G band, for example by ISekii et all 
(l20Q7h . who considered power, phase , and coherence spectra but not spatial 
maps, and bv lNagashima et~al. (|2QQ7h who looked at the power distribution in 
and around a sunspot. In the latter article there was no high-frequency acoustic 
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Figure 2. The mean magnetic inclination (a), from HMI vector magnetograms, courtesy K. 
Hayashi, an AIA 304 A image (b), and an AIA 171 A image (c), for the area of interest on 23 
May 2010. The contours overlaid on panels a and b show the 500 G level of the mean magnetic 
field strength from HMI line-of-sight magnetograms. 



halo detected, but there was a narrow band of excess power surrounding the 
sunspot at all wavelengths in the G band, and pronounced suppression of power 
at both wavelengths in the penumbra and plage, while the umbra showed excess 
po wer in the H line, especially above 4.5 mHz. 

ISchunker and Braunl ( 201lf ) have examined the power distribution for veloc- 
ity observations from MDI and found that the excess high-frequency power 
corresponds to regions with the magnetic-field inclination (as deduced from 
potential-field source-surface extrapolation) in the 40-60 degree range. However, 
the exact origin of these high-frequency haloes remains unclear. 

The relative phase of the Doppler velocity and the intensity in various wave- 
length bands has been studied for more than two decades in both resolved and 
unresolved observations. Early (grou nd-based) spati a lly resolved phase measure- 
ments were made, for example, by iDeubner et Ml (I199QL Il992h : in the latter 



article they used three different wavelengths including the Na D line, formed in 
the lower chromosphere. These measurements revealed a complicat ed pattern of 
phase relationships that was not easily explained by simple models. TMoretti et all 
([20070 made ground-based o bserv ations from the South Pole using the MOTH 



instrument ( Finsterle et all l2004j ). taking high-cadence images in the sodium 



and potassium D lines. The sodium line shows clear evidence of a high-frequency 
halo around an active region. 

In this article we will examine the behaviour and relationships of the different 
observables in a small region of the solar surface containing an active region as 
observed by HMI and AIA, and show how the power and phase of the oscillations 
is affected by the presence of local magnetic fields. For this purpose we will show 
observations of NOAA Active Region 11072 on 23 May 2010. We have also 
examined the data for the same region over the two preceding and following 
days; the results are very similar and for reasons of space are not shown in this 
article. 

In Section [2] we will describe our data and analysis and define the phase and 
coherence spectra. In Section [3] we show the region in each of the observables; 
in Section |4] we present the power maps for each observable, and in Section [5] we 
show the phase and coherence spectra. In Section [6] we discuss our findings. 
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Figure 3. Slices through the smoothed relative-power maps for HMI Iq, at 3 mHz (a), 5 mHz 
(b), 7 mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line 
shown in panels a-d (e). The grey-scale bar applies to all of the panels. 




Figure 4. Slices through the smoothed relative-power maps for HMI V at 3 mHz (a), 5 mHz 
(b), 7 mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line 
shown in panels a-d (e). The grey-scale bar applies to all of the panels. 

2. Data and Analysis 

2.1. Helioseismic Observables from the Solar Dynamics Observatory 

The HMI observables are produced from observations of the Fe I line at 6173 A 
by combining "filtergrams" taken with different configurations of the instru- 
ment's filters. For each 45-second interval there is a line-of-sight magnetogram, 
a Doppler- velocity image [HMI V], a continuum intensity image [HMI Ic], a Line 
Depth image [HMI Ld], and also a Line Width image, which is not considered 
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Figure 5. Slices through the smoothed relative-power maps for HMI at 3 mHz (a), 5 mHz 
(b), 7 mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line 
shown in panels a-d (e). The grey-scale bar applies to all of the panels. 
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Figure 6. Slices through the smoothed relative-power maps for AIA 1700 A, at 3 mHz (a), 5 
mHz (b), 7 mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal 
line shown in panels a-d (e). The grey-scale bar applies to all of the panels. 



here. For this investigation we consider the HMI 1q and HMI V observables and 
also a Line Core intensity observable [HMI II] formed by taking the difference 
between HMI I c and HMI Ld. The pixel size of the 4096 x 4096 pixel images is 
ap proximately - 5", w ith an optical resolution of Y 



Norton et 



(|2006h found that the height of formation of the HMI 6173 A 
line spans the range from 20 km at the wings to 270 km at the core, while 



Fleck, Couvidat, and Strausl ( 201lf ) have calculated that the height of formation 



of the HMI Doppler velocity is approximately 100 km, slightly lower than that of 
the Ni 6768 A line used by MDI. However, the finite resolution of the instruments 
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Figure 7. Slices through the smoothed relative-power maps for AIA 1600 A, at 3 mHz (a), 5 
mHz (b), 7 mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal 
line shown in panels a-d (e). The grey-scale bar applies to all of the panels. 



raises the effective height of formation by about 50 km at disc centre in both 
cases. We can therefore say that the HMI Iq observable is formed considerably 
lower in the atmosphere than HMI with HMI V at an intermediate depth. 
All of these heights apply to quiet Sun, and may well be altered in active regions. 

The AIA 1600 A images for the day analysed have a cadence of 24 seconds, 
while the cadence for 1700 A was 48 seconds. (More recent data have used a 
consistent 24-second cadence for both, but the cadences were adjusted several 
times during the early months of AIA observations.) These images are also taken 
at 4096 x 4096 pixels, but because AIA allows more room around the limb 
of the solar image, the pixel size is approximate ly 0.6" . The 1700 A band is 
believ ed to be formed in the upper photosphere; I Judge, Tarbell and Wilhelm 
(|2001[ ) put it at 300-550 km, around the temperature minimum between the 
upper photosphere and lower chromosphere. The 1600 A band includes the C IV 
line and samples both the transition region at the top of the chromosph ere and 
the upper photosphere. The calculations of Fossum and Carlssonl (|2005f ) for the 
equivalent TRACE bandpasses give mean formation heights of 360 km and 480 
km for the 1700 and 1600 A bands, with FWHM of 385 and 185 km respectively; 
these authors also note that the passbands are double peaked, with the 1700 
A in particular having two almost equal components centred at 402 ± 259 and 
208 ±174 km. The helioseismic signal in the 1600 A band, therefore, is in general 
somewhat contaminated by flare activity, but there were no strong flares during 
th e day analys e d. It seems unlikely that the strong global p-mode signal found 
bv lHowe et all ( 2011 ) in both 1600 A and 1700 A bands could originate as high 
as the transition region, and we therefore suspect that the helioseismic response 
in this bandpass is associated with the continuum rather than the C IV line 
and pertains to a height range not very far above that involved in the 1700 A 
band. We should note that the calibration and registration of the AIA data are 
not optimized for helioseismology in the same way as the HMI data, so there 
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is a possibility of position errors and small drifts, of the order of a fraction of 
a pixel over several days, which could affect comparisons at the pixel level but 
should not be a serious problem for the analysis described here, which involves 
smoothing over several detector pixels. 



The data were obtained from the JSOC web site, http://jsoc.stanford.edu/ 



The website allows the extraction of a "tracked hg_patch" from each image; 
these patches can be chosen to be centred on a given heliographic location but 
are not remapped. For each observable we obtained a 24-hour sequence of such 
patches centred on the location of AR 11072, covering the whole of 23 May 2010. 
These patches were based on images from the aia.levl series and the (level 1.5) 
series hmi.V_45s, hmi.lc_45s, and hmi.l_d_45s; these are the most final calibrations 
of the data available at the time of writing. The data were then remapped to a 
common 201 x 201 pixel grid, 15 heliographic degrees on each side, evenly spaced 
in heliographic latitude and longitude, and interpolated to a uniform 45-second 
cadence. A running mean over 15 minutes was subtracted from each remapped 
pixel time series to remove effects such as daily variations due to the orbit of 
SDO; this will also reduce the reliability of the results at the low- frequency end 
of the spectrum. 

2.2. Phase and Coherence Functions 

The relationship between two time series in the Fourier domain can be expressed 
in terms of the cross-spectrum, the phase difference, and the cross-coherence, 
which we define as follows. 

Let Pa{v) and Pb{v) be the complex Fourier spectra of two time series. The 
"cross spectrum" is defined as 

CROSS(z/) = P A b = P A (v)*Pb(v) (1) 

The phase difference [5(/)a,b] can be written as 

5<t> AB {y) = axg(P AB (v)). (2) 

We adopt the convention that a positive value of 5(/)ab means that A leads B. 
The coherence spectrum, a measure of the correlation of the two spectra as a 
function of frequency, can be written as 

COH AB (,) = - ^BiM . (3) 
V(\Pa{v)\ 2 )(\Pb{v)\ 2 ) 

The angle brackets formally denote the expectation value; in prac tice this can be 
appro ximated by a mean over many short spectra as in the work of lElsworth et a/.^ 



(1994) , a smo o thing by a running mean in the frequenc y domain, as u s ed by 
Jimenez et all (Il999h, or, for resolved- Sun observations ([Severing et all 12001 



Barban, Hill, and KrasL 12004ft an average over all azimuthal orders [m] for a 



given degree [£] . In the present work we use three-dimensional boxcar smoothing 
of the data cube of one-dimensional power spectra, equivalent to taking a running 
mean over three remapped pixels in each spatial direction in turn and then over 
0.3 mHz in temporal frequency. 
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3. Images 

In Figure [1] we show the mean image for each of the observables for 23 May 
2010. On this date the active region was somewhat mature, with a well-defined 
leading spot and a following-polarity region consisting mostly of strong plage 
with a couple of small sunspots, and it was close to the central meridian. In the 
Doppler-velocity panel, the supergranulation is clearly visible; the active region 
disrupts the supergranular pattern and produces the familiar appearance of an 
upflow on one side of the penumbra and a downflow on the other, an artefact of 
the projection of the horizontal Evershed flow on the line of sight. The sunspot 
umbra and penumbra appear as darkened areas in all of the intensity observables, 
but the weaker fields in the following-polarity region and surrounding network 
are bright in the AIA bands and in HMI II . The network fields clearly correspond 
to the supergranule boundaries. For additional context, we show in Figure [2] 
the mean magnetic-field inclination from Milne-Eddington inversions of HMI 
magnetograms and the AIA images at 304 and 171 A, which show the structure 
of the overlying coronal loops that connect the leading and following regions. 



4. Power Maps 

In order to study the local power distribution in each observable, we calculated 
the power spectrum in each pixel of the remapped data cubes. The spectra were 
then smoothed over 25 bins (about 0.3 mHz) in temporal frequency and three 
bins in each spatial direction, and the power at each frequency was normalized 
by the mean power for that frequency across the whole data cube. 

We now examine the power maps for each observable in turn, taking them in 
ascending order of height of formation. 

4.1. HMI I c 

Figure [3] shows four sample relative power maps at 3, 5, 7, and 9 mHz, and 
a slice through the stack of relative-power maps at constant latitude passing 
through the body of the active region, for the HMI Ic data cube. The acoustic 
cut-off frequency in quiet Sun is around 5 mHz. In this case, and for all the other 
intensity observables, we divided each time series by the mean intensity for the 
pixel to remove the unwanted effects of brighter regions. Power is suppressed in 
the active region at all frequencies; there is very little variation in the pattern of 
power suppression with frequency, although there is slightly stronger suppression 
in the five-minute band and above than in the granulation (< 2.5 mHz) range. 
There is no sign of any haloes of excess power. The narrow spike of excess power 
within the sunspot around 5 mHz is seen in only one of the other four days 
examined and may be a noise artifact. The power map shows somewhat finer 
detail at higher frequencies, presumably due to the way that the low- frequency 
limit of the p-mode spectrum at the /-mode ridge shifts towards higher spatial 
frequencies with increasing temporal frequency. 



SOLA: ms_arxiv.tex; 9 August 2012; 0:46; p. 9 



R. Howe et al. 



4.2. HMI Velocity 

In Figure |4] we show the results for the HMI V power. Here we see the well- 
known suppression of the acoustic power in the active region in the five-minute 
frequency band, and a weak, broad halo of excess power surrounding the active 
region at 5 mHz. In the 7 mHz band the excess power has shrunk to a narrow belt 
outlining the active region and surrounded by a wider region of suppressed power, 
with a hint of a more diffuse enhancement beyond that. At this frequency, the 
suppression of power within the sunspot is much weaker and even cancelled out 
by the enhancement around the edges. Examination of the longitude-frequency 
slice suggests that the halo of enhanced power contracts towards the active region 
with increasing frequency. 



4.3. HMI I L 



Figure [5] shows slices through the power maps for HMI II- Here, there is a strong 
halo of excess power - so strong that we have used a wider greyscale range than 
for the Ic and V maps - above about 6 mHz surrounding the active region, but 
farther out than that seen in the velocity observations. This region is beyond 
that where there is any detectable concentration of surface magnetic field, and 
has no obvious correspondence to the morphology of the overlying magnetic 
loops as revealed in the AIA 304 A and 171 A images. At 7 mHz, it occupies 
about the same area as the outer band of reduced power in the velocity power 
map. Within the halo, there are "tendrils" of power suppression extending out 
from the active region that correspond to concentrations of plage field, giving 
it a patchy appearance. As in the velocity observations, this halo appears to 
contract with increasing frequency, and at 9 mHz it has collapsed to a narrow 
band around the sunspots, much the same as that seen in t he HMI V at 7 



mHz. The existence of the halo is consistent with the findings of iJain and Haber 



(|2002[ h who saw a similar phenomenon in the MDI line-depth observations. 



4.4. AIA 1700 and 1600 A 



Finally, we show the power-map slices for AIA 1700 A (Figure [6]) and AIA 
1600 A (Figure [7J). As in HMI II, the power in the active region is more strongly 
suppressed at 5 mHz than in the five-minute band, and suppression is also visible 
for weak network-field concentrations in the quiet Sun. Also, as seen in HMI 
II, there is an outlying halo of high-frequency excess p ower surroundin g the 
active region. This appears to contradict the findings of iMuglachl (2003), who 
did not find any excess power in these bands in TRACE observations. However, 
we reiterate that this power excess is at a greater distance from the active 
region than the halo in velocity, and also occurs at higher temporal frequencies. 
Haloes of excess power have also be en seen in ground-ba sed observations using 
chr omospheric l i nes, fo r example bv lMoretti et all ( 20071 ) in the Na D line and 
byi raun et al. in Ca K. The halo region contracts and weakens with 

increasing frequency, becoming indistinct by about 10 mHz, but the collapse 
and encroachment on the active region seen in the HMI V and HMI II cases do 
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Figure 8. Slices through the HMI J C /HMI V phase, at 3 mHz (a), 5 mHz (b), 7 mHz (c), 9 
mHz (d), and in longitude and temporal frequency along the horizontal line shown in panels 
a-d (e). The colour-scale bar applies to all of the panels. 




Coherence 



Figure 9. Slices through the HMI I c /HMI V coherence, at 3 mHz (a), 5 mHz (b), 7 mHz 
(c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line shown in 
panels a-d (e). The grey-scale bar applies to all of the panels. 



not happen up to the 10.4 mHz Nyquist frequency of the 1700 A data set, and 
a re-analysis of the 1600 A data with the full 24-second temporal resolution did 
not reveal any such effect at higher frequencies. 



5. Two-Dimensional Phase and Coherence maps 

The phase and coherence of the individual remapped pixel time series, as a 
function of position and frequency, can be derived using Equations (|2|) and (j3]). 
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Figure 10. Slices through the HMI I L /RMl V phase, at 3 mHz (a), 5 mHz (b), 7 mHz (c), 9 
mHz (d), and in longitude and temporal frequency along the horizontal line shown in panels 
a-d (e). The colour-scale bar applies to all of the panels. 
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Figure 11. Slices through the HMI I L /YLMl V coherence at 3 mHz (a), 5 mHz (b), 7 mHz 
(c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line shown in 
panels a-d (e). The grey-scale bar applies to all of the panels. 



When the coherence is low there is little relationship between the two signals 
and we can expect the apparent phase difference to be close to zero after spatial 
averaging. 

In Figures [8] to [15] we show these quantities for the combination of HMI V 
with HMI I c , HMI J L , AIA 1700 A, and AIA 1600 A. The velocity is treated 
as positive in the direction away from the observer. In all cases, the coherence 
with the velocity is sharply reduced within the active region, which is under- 
standable due to the greatly reduced mode power and consequent reduction 
in signal-to- noise ratio. There is, however, a local increase in coherence in the 
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Figure 12. Slices through the AIA 1700 A/HMI V phase at 3 mHz (a), 5 mHz (b), 7 mHz 
(c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line shown in 
panels a-d (e). The colour-scale bar applies to all of the panels. 
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Figure 13. Slices through the AIA 1700 A/HMI V coherence at 3 mHz (a), 5 mHz (b), 7 
mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line shown 
in panels a-d (e). The grey-scale bar applies to all of the panels. 



core of the sunspot for all three variables. This may be a real effect, but we 
should be cautious of common spurious, or at least complicating, effects within 
the sunspot umbra, such as noise or tracking errors interacting with the steep 
gradients of power and intensity, contamination by other spectral lines due to 
reduced temperatures, and possible cross-talk between the HMI observables in 
the presence of strong fields. 

The phase difference between HMI V and HMI Iq shows a sharp discontinuity 
at around 2.5 mHz, the upper limit of the granulati on signal; thi s phenomenon 
has been well known since the work of, for example, IStaigerl (|1985f ) as illustrated 
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Figure 14. Slices through the AIA 1600 A/HMI V phase at 3 mHz (a), 5 mHz (b), 7 mHz 
(c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line shown in 
panels a-d (e). The colour-scale bar applies to all of the panels. 
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Figure 15. Slices through the AIA 1600 A/HMI V coherence at 3 mHz (a), 5 mHz (b), 7 
mHz (c), 9 mHz (d), and in longitude and temporal frequency along the horizontal line shown 
in panels a-d (e). The grey-scale bar applies to all of the panels. 



bv lDeubner et all ( 1990f ), and of I Jimenez et all ( 19991 ). In the core of the active 
region this boundary is shifted to lower frequencies, perhaps due to the sup- 
pression of convective motions by the surface fields. The weaker fields of the 
plage area have very little effect on the Ic~V phase in the five-minute band, 
but do modulate the Il~V phase difference. In the HMI II and AIA cases, the 
discontinuity is smaller and the low-frequency behaviour more complicated, with 
the phase in the plage region heavily disrupted. 
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While the quiet-Sun phase difference between HMI 1q and HMI V becomes 
slowly less negative with increasing frequency, the opposite is true for HMI II 
and the UV bands. 

Another interesting feature is that the high-frequency coherence between the 
HMI V and the UV intensities is higher in the "halo" zone around the active 
region than in the surrounding quiet Sun; this region (which, as mentioned above, 
does not correspond to any detectable magnetic fields) also shows a shift in the 
phase between these observables that varies smoothly with distance from the 
active region, going from a difference of about —100 degrees for the 1700 A and 
— 160 degrees for the 1600 A band in quiet Sun far away from the active region 
at 7 mHz to nearly in-phase within it. While the coherence of HMI V with 
AIA 1600 A and 1700 A look very similar at lower frequencies, the 9-mHz map 
shows a different pattern; while for the 1700 A band the coherence is weakly 
enhanced over most of the extended halo region except for the active region 
core and an external fringe of reduced coherence, for the 1600 A band there is a 
high-coherence fringe surrounding a patch of reduced coherence over the wider 
active region. 

On the other hand, the coherence between HMI V and HMI II is reduced in 
the same region and frequency range, and the phase difference varies from —130 
degrees in quiet Sun to —90 degrees in the active region. At 3 mHz, conversely, 
the phase difference between HMI V and both HMI I L and AIA 1600 A and 
1700 A is more negative within the active region than in quiet Sun. It is curious 
that the quiet- Sun phase is more negative than that for the AIA bands, as this 
would seem to imply (if we were dealing with upward-traveling waves) that the 
HMI II is actually probing a higher level of the atmosphere than the UV bands. 

6. Discussion 

We have examined the spatio-temporal power distribution around an active 
region in a number of HMI and AIA observables, and the phase and coherence 
relationships between the intensity observables and the HMI Doppler velocity. 

Five-minute power in all observables is suppressed in the sunspot (which is 
dark at all wavelengths) and also in plage areas (bright in AIA bands and HMI 
II)- Above the acoustic cut-off frequency the behaviour is more complicated. 
Power in HMI Iq is suppressed in the presence of surface fields at all frequencies, 
while power in HMI V is enhanced in a narrow zone around field concentrations 
(especially plage) and suppressed in a wider surrounding area; these results are 
consistent with earlier work using MDI. Power in HMI II and the AIA bands is 
suppr essed in areas of surface field but - in contrast to the results of iMuglachl 
(|2003h - enhanced in an extended area around the active region. In the HMI 
II case, the pattern of narrow halo, encroaching on the region where the five- 
minute power is suppressed, and surrounding region of suppressed power, seen 
for the HMI V at 7 mHz, is seen instead at 9 mHz. For the UV bands, however, 
this pattern does not appear at any frequency; instead, the halo power fades 
with increasing frequency and vanishes above about 10 mHz while the power 
suppression remains. 
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In all cases but HMI /c, the regions of enhancement and suppression of power 
appear to move inwards towards the active region at increasing frequency. 

The relative phase of the observables is altered around active regions. While 
the apparently in-phase behaviour of most observables in strong-field areas is 
associated with low coherence and thus of low significance, there are exceptions 
to this. In particular, in the 7 mHz frequency band there are areas close to the 
active region where the AIA 1700 A and 1600 A bands are both close in phase to 
the HMI velocity and coherent with it, whereas in the quiet Sun well away from 
the active region the phase difference is close to 90 degrees and the coherence is 
lower. This effect is not seen in the HMI II, which is otherwise very similar in 
its behaviour to the UV bands, although there is a phase difference between the 
two observables. At 9 mHz, the 1600 and 1700 A bands also show qualitatively 
different behaviour, suggesting that though there may be substantial overlap 
between the two wavelength bands, the highest-frequency waves are sensing 
different layers of the atmosphere. 

Clearly, within the "halo" zone surrounding the active region the propa- 
gation and reflection of the high-frequency waves differs from that over quiet 
Sun, whether due to alteration in the height of formation of the observables, 
thermal changes, or the direct effect of magnetic fields in trapping, scattering 
or transforming the waves; also, there are qualitative differences between the 
purely photospheric observables - even at the relatively high formation level of 
the HMI I L - and the UV bands. 

In general, the higher the observable is formed, the higher the frequency at 
which the halo of excess power is found, and the more extensive the halo at 
the lowest frequency at which it appears. It is tempting to interpret these halo 
effects in terms of spreading magn etic-field canopies, as has been suggested, 
for example, bv lMoretti et all ( 2007 ), but we note that the extent of the high- 
frequency halo does not exactly match the morphology of the overlying field 
seen in the EUV images. The size of the halo and the inner and outer regions 
of suppressed power appear to vary with frequency as well as with the height of 
formation of the observable used, contracting with increasing frequency. If the 
patterns of power distribution are related to spreading out of magnetic fields 
with height, this would imply that the higher- frequency pseudomodes are being 
reflected or absorbed at lower layers than the lower- frequency ones. The heights 
of formation themselves may also be altered in the presence of magnetic fields, 
further complicating the picture. 

These effects no doubt also influence the changes in the local helioseismic 
parameters of modes ab ove the acoustic cut-off frequency, seen for example by 
Howe et~al\ ( 20041 . 2008 ). which tend to be in the opposite sense from those 



experienced by modes in the five-minute range. 

Eventually, we hope to carry out multi- wavelength helioseismic analysis using 
HMI and AIA data together to extract parameters for local and global acoustic 
spectra. As the results described here demonstrate, any such analysis will need 
to take into account the effects of active regions on the phase and coherence of 
the oscillations. 

The continuous, high-cadence, full-disc observations of HMI and AIA allow us 
to study the behaviour of waves in the photosphere and lower chromosphere at 
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a level of detail that has not previously been possible. These observations may 
help us to improve our understanding of the interaction of waves and magnetic 
fields in the different layers of the photosphere. While the analysis described 
here does not distinguish between true waves and other short-period spatio- 
temporal variations, we hope in the future to examine the three-dimensional 
power spectrum and its phase and coherence in order to separate the waves 
from granulation and other background effects. 
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